AD/A-006  429 


A GENERALIZED  CAP  MODEL  FOR  GEOLOGICAL 
MATERIALS 


I.  S.  Sandler,  et  al 


Weidlinger  Associates 


Prepared  for: 

Defense  Nuclear  Agency 


14  November  1574 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


S' 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  page  CHJ,rn  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

J.  ' REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

DNA  3443T 

3.  RECIPIENT'S  CATALOG  NUMBER 

ADjA  -op  6 54?? 

4.  TITLE  (and  Subtitle) 

A GENERALIZED  CAP  MODEL  FOR  GEOLOGICAL 
MATERIALS 

5.  TYP^Or  REPORT  ft  PERIOD  COVERED 

Topical  Report 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORfsJ 

I.  S.  Sandler 

F.  L.  DiMaggio 

G.  Y.  Baladi 

8.  CONTRACT  OR  GRANT  NUMBERfsJ 

DNA  001-73- C-0023 

9.  PERFORMING  ORGANIZATION  NAME  AND  AOORE5S 

Weidlinger  Associates 
110  East  59th  Street 
New  York,  New  York  10022 

!0.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  a WORK  UNIT  NUMBERS 

NWED  Subtask 
Y99QAXSB047-66 

II.  CONTROLLING  OFFICE  NAME  AND  AODRESS 

Director 

Defense  Nuclear  Agency 
Washington,  D.  C.  20305 

12.  report  date 
14  November  1974 

13.  NUMBER  OF  PAGES 

■if 

14  MONITORING  AGENCY  NAME  A ADDRESSfJ/  dlllerent  from  Controlling  O flier) 

IS.  SECURITY  CLASS  (of  this  report) 

UNCLASSIFIED 

15  a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

•6.  DISTRIBUTION  STATEMENT  ( of  this  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  ebstrnct  entered  in  Block  20,  II  dillcrcnt  from  Report) 

IB.  SUPPLEMENTARY  NOTES 

This  work  sponsored  by  the  Defense  Nuclear  Agency  under  Subtask 
Y99QAXSB047-66. 

19.  KEY  WORDS  /Conf/mje  on  reverse  side  if  necessary  and  identify  by  block  number) 

Geology  Dilatancy 

Cap  Model 
\ Plasticity 

| Stress-Strain  Curves  PRICES  SUBJECT  10  CHANGE 

20.  ABSTRACT  (Continue  on  reverse  side  if  necessary  and  identify  by  block  numbet) 

; The  full  scope  and  theoretical  basis  of  the  cap  model  is  described  so  as  to 

, demonstrate  and  define  the  capabilities  and  limitations  of  the  model.  A cap  model 
falls  within  the  framework  of  the  classical  incremental  theory  of  plasticity  and  is 
based  on  a loading  function  which  serves  as  both  a yield  surface  and  plastic  potential. 
The  model  describes  material  behavior  in  compression  and,  in  general,  some  type 
of  tension  behavior  must  be  postulated  for  completeness.  Tensile  behavior  will  not 

DD  I jiuw  1473  £01’  Reproduced  by  / TITvIPT  ASST  VTF  D 

J NATIONAL  TECHNICAL  / UNlLAbbll  1LU 

INFORMATION  SERVICE  ecurity  classification  of  this  page  (Wncn  Data 


US  Department  of  Commerc* 
Springfield,  VA.  22151 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  THIS  PAGE(H?imi  Dor*  Entered) 


20:  ABSTRACT  (Continued). 

be  considered  here.  Particular  forms  of  the  cap  model  are  adequate  for  many 
purposes,  however,  it  is  desirable  to  describe  the  model  in  its  most  general  form 
. so  as  to  clearly  indicate  the  adaptability  and  flexibility  as  well  as  the  limits  of  the 
cap  model  approach.  This  has  been  done,  and  the  procedure  for  fitting  the  cap  model 
is  briefly  described. 


" f * t 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfHVicn  Data  Entered) 


TABLE  OF  CONTENTS 


Page 


I INTRODUCTION . 3 


II  GENERALIZED  CAP  MODEL G 


III  PROCEDURE  FOR  FITTING  OF  CAP  MODEL 11 


IV.  AN  EXAMPLE  OF  A GENERALIZED  CAP  MODEL 14 

REFERENCES 22 

APPENDIX  A 23 

APPENDIX  B 25 


1 


.1 


INTRODUCTION 


In  two  previous  studies.  Ref.  [1]  and  [2],  particular 
cap  models  for  soils  and  rocks  were  introduced.  The  purpose 
of  this  report  is  to  indicate  the  manner  in  which  previous 
* models  may  be  generalized,  if  necessary,  to  more  adequately 

describe  the  behavior  of  geological  materials.  The  full 
scope  and  theoretical  basis  of  the  cap  model  is  described  so 
as  to  demonstrate  and  define  the  capabilities  and  limits  of 
the  model.  Also,  the  fitting  procedure  and  use  of  a 
particular  version  of  the  general  model  is  illustrated  by 
means  of  an  example. 

From  a general  point  of  view,  a cap  model  falls  within 
the  framework  of  the  classical  incremental  theory  of 
. plasticity  and  is  based  on  a loading  function  which  serves  as 

both  a yield  surface  and  plastic  potential.  Typically,  the 

'm 

loading  function  is  assumed  to  be  isotropic  and  to  consist  of 
two  parts:  a modified  Drucker-Prager , Ref.  [3],  yield 

condition,  denoted  by 

fl<Jl  , J^)  = 0 (1) 

in  which  and  J 2 are  the  first  and  second  invariants  of  the 
stress  and  deviatoric  stress  tensors,  respectively , together  with 
hardening  plastic  cap 

f2(Jx  , , K)  = 0 (2) 

4 

which  may  expand  or  contract  as  the  hardening  parameter  K 

1 

increases  or  decreases.  These  are  illustrated  in  Fig.  1. 

The  model  describes  material  behavior  in  compression  (J^ £ 0) 
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and,  in  general,  some  type  of  tension  behavior  must  be 
postulated  for  completeness.  Tensile  behavior  (for  which 
J.,  > 0)  will  not  be  considered  here. 

In  the  previous  models  for  soils,  the  hardening 
parameter  was  taken  to  be 


K = 


(3) 


in  which  is  the  plastic  volumetric  strain,  while  for 
rocks 


,t. 


K ~ 


' J2> 


• p ^ . D 2 • D 2 

(ep  + (ep  + (ep 


d t 


(M 


.p 

in  which  the  are  the  principal  components  of  the  plastic 
strain  rate  tensor  and  t is  time.  The  use  of  Eq.  (3) 
permitted  the  cap  to  reverse  itself  when  a point  on  the  yield 
curve  = 0 was  reached,  thus  controlling  the  excessive 
dilatancy  predicted  for  soils  by  the  Frager-Drucker  model. 

The  use  of  Eq.  (4)  for  rocks,  which  does  not  permit  the  cap 
to  move  back,  ensures  dilato.ncy  v.hile  permitting 
hysteresis  in  a hydrostatic  load  - unload  cycle.  Such 
volumetric  hysteresis  is  not  present  in  the  Prager-Drucker 
model . 

In  the  earlier  models,  che  elastic  portion  of  the 
behavior  (which  is  most  important  in  determining  the  unloading 
and  reloading  behavior  of  the  model)  was  generally  assumed 
to  be  linear,  i.e.,  it  was  described  by  constant  bulk  and 
shear  moduli. 


r 


f 

I 
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It  was  found  that  good  fits  (well  within  the  scatter  of 
the  experimental  data)  of  stress-strain  curves  and  loading 
paths  were  obtained  using  an  exponential  function  for  f^  and 
an  ellipse  (tangent  to  f^  = 0 for  rocks,,  and  with  a horizontal 
tangent  at  its  intersection  with  f^  = 0 for  soils)  for  f£  = 0. 

While  these  particular  forms  of  cap  model  are  adequate 
for  many  purposes,  it  is  desirable  to  describe  the  model  in 
its  most  general  form  so  as  to  clearly  indicate  the  adapt- 
ability and  flexibility  as  well  as  the  limits  of  the  cap 
model  approach.  This  is  done  in  Section  II.  In  Section  III, 
the  procedure  for  fitting  the  cap  model  is  briefly  described, 
while  Section  IV  gives  an  example  of  a particular  form  of  cap 
model  which  was  used  in  ground  shock  calculations. 
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II  THE  GENERALIZED  CaP  MODEL 


The  classical  theory  of  plasticity  allows  for  a broad 
range  of  material  behavior,  and  the  cap  model  falls  within 
this  range.  Many  previous  applications  of  plasticity  theory 
for  metals  involved  the  assumption  that  volumetric  strains 
are  purely  elastic.  The  cap  model  however,  is  predicated 
on  the  fact  that  the  volumetric  hysteresis  exhibited  by  many 
geologic  materials  can  also  be  described  by  a plasticity 
model,  if  the  model  is  based  on  a hardening  yield  surface 
which  includes  conditions  of  hydrostatic  stress.  Guidelines 
as  _o  how  this  may  be  done  have  been  provided  by  Drucker, 

Ref.  [4],  whose  stability  postulate  is  sufficient,  although 
not  necessary,  to  satisfy  all  thermodynamic  and  continuity 
requirements  for  continuum  models.  Stability  ensures  that  all 
physically  reasonable  initial-boundary  value  problems  are  pro- 
perly posed  in  the  mathematical  sense. 

The  basic  implications  of  Ref.  [4]  with  respect  to  the 
plastic  portion  of  the  material  behavior  are: 


1)  The  yield  surface  should  be  convex  in  stress  space 

2)  The  loading  function  and  plastic  potential  should 
coincide  (associated  flow  rule) 

3)  Plastic  strain  or  work  "softening"  should  not  occur. 


in  which  a,,  and  . are  the  components  of  the  stress 

J * J 


and  plastic  strain  tensors. 

These  conditions  allow  considerable  leeway  in  choosing  the 
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functional  forms  f^  and  f^  in  Eqs.  (1,  2),  but  the  particular 
functions  given  in  Section  -IV  seem  capable  of  adequately 
describing  most  geological  materials. 

It  should  be  noted  that  work  by  Mroz,  Ref.  [5]  and  Bleich, 
Ref.  [6],  suggest  that  it  may  be  possible  to  violate  the  above 
conditions  in  certain  cases  without  destroying  the  stability 
of  the  model. 

The  fitting,  of  stress-strain  curves  and  loading  paths 
using  a cap  model  can  be  mproved  by  introducing  a nonlinear 
elastic  component  of  behavior.  This  has  been  done  by  replacing 
the  constant  bulk  and  shear  moduli  in  the  linear  elastic 
stress  increment-strain  increment  equations 


dJ  = 3K  dl 

(6) 

ds.  . = 2G  de. . 

(7) 

by 

II 

c-< 
>— * 

(8) 

and 

G = GU^) 

(9) 

In  the  above  equations  I is  the  first  invariant  of 

the  strain 

tensor  and  s..  and  e..  are  the  stress  and  strain  deviators, 
il  il 

respectively . 

Additional  flexibility  in  fitting  experimental 

data  can 

be  introduced  into  the  model  by  generalizing  Eqs.  (8)  and 

(9)  to 

K = K ( J L , k) 

(10) 

0 = G ( J ^ > k) 

(ID 
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It  should  be  noted  that  Eqs.  (10)  and  (11),  or  Eqs.  (8) 
and  (9),  correspond  to  a hypoelastic  model.  Ref.  [7],  with 
a positive  definite  elastic  internal  energy,  W,  which  is 
independent  of  stress  path.  This  may  be  shown  by  writing 


(C  . . 
13 


W = 


(O.  . 

ij 


a . . de  . . = 
13  13 


j 1 ^ s ..  ^ 3 1 

(sii  + 3 J1  6ii)(  2G  + TF  fiij* 


s . . 

13  1 

■7TZ  s . .ds  . . + 
2G  ij  13 


J2  dJ^ 


W JidJi 


1 d(J1") 


(12) 


2G(J2  , K)  J 18K(J1  , K) 


During  elastic  deformation  K is  constant,  so  that  the  integrals 
in  the  last  member  of  Eq.  (12)  depend  only  on  the  current 

t 

values  of  J2  and  as  well  as  on  K . Therefore,  W is 
independent  of  path  during  purely  elastic  deformations. 

t 

Further,  since  G and  K are  always  positive,  as  is  J2  , W is 
positive  definite.  Therefore,  there  is  no  possibility  of 
energy  generation  by  the  model. 


In  problems  involving  wave  propagation  at  sites  consis- 
ting of  layers  of  both  soil  and  rock,  it  is  desirable  to  be 
able  to  use  the  same  digital  computer  program  for  layers  of 
both  kinds.  This  is  not  possible  when  the  different  harden- 
ing parameters  of  Eqs.  (3)  and  (4)  are  used.  If,  instead  of 
Eq.  (4),  the  hardening  parameter  for  rocks  is  taken  to  be 


> 


8 


the  maximum  previous  value  of  plastic  compaction 


tc  = ep.  . ( 

v (max) 

a generalized  program  which  may  be  used  for  both  soils  and 
rocks  can  be  easily  written.  It  should  be  pointed  out  that 
the  use  of  Eq.  (13)  instead  of  Eq.  (4)  does  not  permit  the 
introduction  of  dilatancy  until  the  stress  point  first 
reaches  the  yield  surface  f^  = 0,  but  this  loss  in  accuracy 
of  the  fit  may  be  less  important  than  the  improved  efficiency 
of  using  one  generalized  model  for  all  geological  materials. 

It  is  also  noteworthy  that  in  those  problems  in  which 
no  advantage  is  gained  by  using  essentially  the  same  compu- 
tar  code  for  rock  as  that  for  soil,  and  where  the  rock  is 
not  too  porous,  a model  incorporating  a yield  surface  f^  = 0, 
no  cap,  and  different  bulk  moduli  for  loading  and  unloading, 
is  sufficiently  accurate  and  simpler  to  use  than  any  cap 
model.  This  possibility  is  discussed  in  more  detail  in 
Appendix  A. 

It  is  also  possible  to  include  other  features  in  cap 
models,  such  as  anistropy,  rate  dependence,  and  hardening 
of  the  modified  Drucker-Prager  portion  of  the  yield  surface 
(isotropic  and/or  kinematic  hardening).  Further,  in  some 
cases  it  may  be  desirable  to  replace  hypoelastic  behavior 
with  hyperelastic  behavior.  Ref.  [7],  within  the  yield  surface. 
Since  these  additional  possibilities  have  not  been  sufficiently 
studied,  they  will  not  be  included  in  this  paper,  but  will 


be  reported  upon  later. 

If  it  is  necessary  to  introduce  finite  strains,  Jaumann's 
stress  rate  and  the  rate  of  deformation  replace  the  stress 
and  strain  rates  of  Eqs.  (6)  and  (7),  Ref.  [7]. 
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Ill  PROCEDURE  FOR  FITTING  OF  THE  CAP  MODEL 


The  procedure  for  obtaining  the  functional  forms  and 
parameters  used  in  cap  models  constructed  for  use  in  ground 
shock  computations  is  based  on  representative  material  data 
obtained  from  laboratory  tests  on  material  samples. 

Recently,  in-situ  material  tests  have  come  to  be  recognized 
as  having  important  potential  in  ascertaining  material 
behavior  in  situations  of  importance  in  ground  shod, 
computations.  Work  is  currently  in  progress  along  several 
fronts  to  determine  methods  by  vhich  such  tests  can  be 
incorporated  into  the  overall  procedure  for  determining 
material  behavior. 

The  representative  material  behavior  generally  consists 
of  uniaxial  strain  and  triaxial  stress  data.  Sometimes, 
hydrostatic,  proportional  loading,  direct  shear  and/or  other 
tests  are  available.  The  first  step  in  the  fitting  procedure 
is  to  employ  the  unloading  portion  of  these  cests  to 
determ.ne  appropriate  elastic  behavior  of  the  model,  since  the 
cap  model  behaves  elastically  during  initial  unloading  in 
these  tests.  For  example,  as  long  as  the  model  behaves 
elastically,  unloading  behavior  indicates  the  bulk  modulus  K 
in  hydrostatic  tests,  the  shear  modulus  G in  triaxial  stress 
tests,  and  the  combination  K + 4/3  G in  uniaxial  stress  tests. 
Other  tests,  if  available,  may  be  used  to  check  or  adjust  the 
overall  fit. 

The  next  step  in  the  fitting  procedure  is  to  establish 
the  failure  envelope,  i.e.,  the  portion  of  the  yield  surface 
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which  limits  the  shearing  stresses  that  the  material  can 
withstand.  While  the  failure  envelope  could  be  chosen  as  a 
work  or  strain  hardening  yield  surface,  it  is  generally 
adequate  and  much  simpler  to  assume  it  to  be  ideally  plastic. 
The  failure  envelope  is  generally  obtained  using  failure  data 
from  triaxial  stress  and  proportional  loading  tests.  This 
data  is  fit  by  a function  of  the  stresses,  and  is  usually 
assumed  to  involve  only  the  first  stress  invariant  and  the 
second  invariant  of  the  stress  deviators. 

The  remaining  step  in  the  fitting  procedure  is  the  most 
difficult.  The  cap  portion  of  the  model  is  obtained  by  a 
trial  and  error  procedure  in  which  a cap  shape  and  hardening 
rule  are  assumed  and  the  behavior  of  this  assumed  model  is 
computed  and  compared  to  the  representative  material  data. 

If  the  fit  requires  improvement,  a new  set  of  parameters  is 
tried  and  the  procedure  is  repeated.  The  computation  of  the 
model  behavior  can  be  based  on  the  equations  describing  the 
relations  between  the  stress  and  strain  increments  during 
the  common  laboratory  loading  paths.  These  equations  are 
derived  in  Appendix  B for  the  case  of  uniaxial  strain. 

Obviously,  the  success  of  such  a trial  and  error 
procedure  and  the  rapidity  with  which  it  converges  is 
strongly  dependent  on  the  experience  of  the  modeler. 

Knowledge  of  the  effect  on  the  model  behavior  of  changes  in 
the  cap  model  parameters  is  important  for  rapidly  obtaining  a 
satisfactory  fit.  For  example,  the  fitting  procedure  is 
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greatly  simplified  by  the  knowledge  (obtained  through 
experience  in  fitting  several  cap  models  by  trial  and  error) 
that  the  hardening  rule  strongly  affects  the  stress-strain 
curves  for  uniaxial  strain  and  hydrostatic  paths,  while  the 
shape  of  the  cap  plays  an  important  role  in  determining  the 
stress-strain  behavior  for  triaxial  stress  situations  and 
the  stress  path  for  uniaxial  strain.  In  fact,  the  hardening 
rule  has  been  obtained  for  the  most  recent  cap  models  by  using 
a separate  program  to  compute  the  plastic  volumetric  strain 
during  hydrostatic  loading. 


IV  AN  EXAMPLE  OF  A GENERALIZED  CAP. MODEL 


In  this  section  an  example  of  the  type  of  cap  model  used 
in  recent  ground  shock  calculations  is  presented.  Specific 
model  parameters  and  the  resulting  model  behavior  as  well  as 
the  original  material  properties  to  which  the  model  was 
fitted  are  also  presented. 


V 


F(Jl‘VJ2’ev)  = < 


The  failure  envelope  and  cap,  Fig.  1,  which  describe 
the  entire  yield  surface  of  the  material,  intersect  at  = L 
The  complete  yield  function,  F(J,  , 

Fi(ji,yr>  =v^2~~ [A_ce  i]  = ° if  l<ji 

[ F2(J1’*\/^2’£v)  = (VL>2+r2V(X-D2 

= 0 if  L>Ji 

in  which  A,  B and  C are  material  constants.  The  two  quantities 
L and  X are  functions  of  and  represent  respectively 
the  values  of  at  the  center  of  the  cap  and  at  the  intersection 
of  the  cap  with  the  J^-axis.  These  two  quantities  are 
related  to  each  other  through  the  parameter  Z by 

B £ 

Z - X = R [A  - Ce  3 

U if  % < 0 


(14) 


(15) 


L = 


(16) 


0 if  Z > 0 


and  are  related  to  the  plastic  volumetric  strain  by  means 


of  the  hardening  function 
£P  = w[eDX  _ i - 


(D,X-D,f)  2 D X 

aDXe  x * ] -WpX  e 


(17) 
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in  which  W,  D,  a,  D.  , D„  , W„  and  D_  are  material  constants. 
The  quantity  R in  Eqs.  (15,  16)  represents  the  ratio  of  the 
major  to  the  minor  axis  of  the  cap  and  is  given  by 
R R L -R  (L+R  )2 

R(L)  = T&l  [1  + V 1 + R3e 

in  which  R , R , R , R_, R.  and  R_  are  material  constants. 

U 1 / j 4 5 

The  elastic  portion  of  the  model  behavior  is  represented 
by  the  bulk  and  shear  moduli 


K(Jr  £y)  = K 


P ^ — K 1 + ^ 
El 


K1J1 


K 


G(  J2.ej)  ~ gei 


1 + Be^1  2 cosh<K2£v> 

yfiL  Wv 


/ — • p 

-G.-\  / J,  2 cosh  (G  ev) 

1 + ne  V 


in  which  K^,  8,  6,  K^,  K2 , Kg,  y,  G£I , n,  G^  G2>  G 3 and 

G are  material  constants, 
s 


For  the  material  constants  listed  below  the  model  behavior 
is  shown  in  Figs.  2-4  together  with  the  material  representa- 
tive properties  to  which  the  model  was  fitted.  In  general,  the 
model  agrees  quite  well  with  the  representative  properties. 

The  model  parameters  are 


(18) 


(19) 


(20) 
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AXIAL  STRAIN  (%) 


STRESS- STRAIN  BEHAVIOR  OF  MODEL 
VS.  REPRESENTATIVE  MATERIAL 
IN  UNIAXIAL  STRAIN 


REPRESENTATIVE  BEHAVIOR 
MODEL  BEHAVIOR 


FIG.  3 STRESS  STRAIN  BEHAVIOR  OF  MODEL  VS.  REPRESENTATIVE  MATERIAL 

IN  UNIAXIAL  STRAIN 
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STRESS  PATH  BEHAVIOR  OF  MODEL  VS.  REPRESENTATIVE  MATERIAL 

IN  UNIAXIAL  STRAIN 


A = 0.405  ksi 
B = 0.5  ksi-1 
C = 0.403  ksi 
W = 0.101 
D = 2.0  ksi"1 
a = 1.0 

= 16.0  ksi  1 
D2  = 0.0  ksi-2 
= -0.43  ksi-1 
Df  = 3.65  ks"  1 

Ro  = 4-3 

R.l  = 0,8 

R2  = 4.0  ksi  1 

R3  = -3.5 

y = -0.7 


8 
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G 

s 


n 


0.03  ksi 

7.0  ksi 
444.  ksi 
26.6  ksi 
5.28 
-0.924 
0.693  ksi 
80 

267  ksi 
5.4  ksi 
4 

0.0  ksi-1 

3.0  ksi-1 
200 


It  should  be  noted  that, in  general,  the  choice  of  material 
parameters  cannot  be  made  in  completely  arbitrary  fashion. 

For  example, 


A > C 
B 0 
R > 0 


(22) 

(23) 

(24) 
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so  the  parameters  Rg,  R^»  R2»  R^»  and  R^  must  be  chosen 

so  as  to  satisfy  Eq.  (24)  for  all  possible  states  of  the 

material.  Further,  the  loading  condition  in  conjunction 

with  the  yield  condition  F2  = C requries 

3F„  3F 

— - dep  = - tt-~  da. . < 0 

3ep  v 3°ij  13  " 

v J 


and  K > 0 and  G > 0 must  also  be  satisfied  for  all  achiev- 
able stress  states  of  the  material. 
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APPENDIX  A 

It  has  been  found  while  making  fits  for  rocks  of 
low  porosity,  that  the  cap,  = 0,  may  be  chosen  as 
essentially  vertical  except  near  its  intersection  with 
the  yield  surface  f^  = 0.  The  behavior  of  a cap  model 
with  a vortical  cap  and  a bulk  modulus  K (which  may  b.*  a 
constant  or  function  of  pressure)  which  is  the  same  for 
loading  and  unloading  (see  Fig.  5),  i.e.. 


KL  = Ku 


(A-l) 


II 


is  identical  to  the  uncapped  model  (see  Fig.  6)  with 


KL  < Ku 


That  this  is  so  can  be  readily  seen  by  noticing  that  if  an 
associated  flow  rule  is  applied  to  a vertical  cap  only 
plastic  volume  changes  occur.  The  resulting  hysteresis 
in  the  stress- strain  curves  is  also  produced  without  the 
cap  if  Eq.  (A-2)  is  used. 

It  should  be  noted  that  use  of  an  inequality  like  Eq . 
(A-2)  involving  the  shear  modulus  would  violate  the  con- 
tinuity requirements  that  load  paths  infinitesimally  close 
to  a neutral  loading  path  should  result  in  essentially  the 
same  stress-strain  curves. 


(A-2) 


Xi  ' 


23 


-J, 

FIG. 6 IDEALLY  PLASTIC  MODEL  (WITHOUT  CAP) 

WITH  DIFFERENT  BULK  MODULI  ON  LOADING 
AND  UNLOADING 


APPENDIX  B 


CAP  MODEL  BEHAVIOR  IN  UNIAXIAL  STRAIN 


The  general  material  model  considered  here  can  be  described 
in  three  dimensional  Cartesian  coordinates  by  the  incremental 
relations 


d£. . — d£ . . + d£?. 
13  13  13 


(B-l) 


1 - 


1 . 


= or  + 2G  ” -J 


ij 
F (c 


9K  ij 


13 


13 


'kkJ 


13 


’ 

d £?.  * 

*3  \ 


1° 


■v 

< 0 

dA 

3F 

if  F 

30 

ij 

0 

if  F 

(B-2) 


(B-3) 


(B-4) 


in  which  the  summation  convention  has  been  adopted,  5 is 

E P ij 

the  Kronecker  delta,  and  the  d£..  , d£..  and  d£..  denote  the 

13  ij  13 

increments  of  total  strain,  elastic  strain,  and  plastic 
strain,  respectively.  The  do_  denote  the  stress  increments, 
and  dA  is  a coefficient  which  is  non-zero  only  when  plastic 
deformations  occur.  During  plastic  deformation  Eqs.  (B-3,  4) 
become 


F(o.  . , e-  . ) = 0 
13  iJ 


. P 3F 

de.  . = dA  -r— — 
ij  do . . 

J 13 


(B-5) 

(B-6) 


The  elastic  and  plastic  strains  can  be  eliminated  from  Eqs. 
(B-l,  2,  3,  4)  by  differentiating  Eq.  (B-5) 


0 


( B—  7 ) 
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and  substituting  Eq.  (B-2,  6)  into  Eqs.  (B-l,  7)  to  obtain 


d£.  . 

13 


_Lx  , 

6G  ij 


da 


kk 


do . . + dX 
13 


3F 


3a. . 

13 


3F 


3o.  . 

13 


da. . + dX 

13 


3F 


3a. . 

13 


3F 

3ep 

ij 


0 


Multiplication  of  Eq . (B-8)  by  5 gives 


dOj^  3k(d£j,^  — dX 


kk 


3a  rs 
rs 


(B-8) 

(B-9) 


(B-10) 


and  substitution  of  Eq.  (B-10)  into  Eq.  (B-8)  leads  to 


da..  = 2Gde . . + (K  - —■)&.  . (d£..  - dX-r^-  5 ) 

13  xj  3 13  kk  3a  rs 


rs 


- 2G  dX 


_3_F_ 

3a.  . 

13 


Multiplication  of  Eq.  (B-ll)  by 


3F 

3a.  . 

3-3 


and  subtraction  of 


the  resulting  relation  from  Eq.  (B-9)  gives  an  equation 
which  may  be  solved  for  dX 


dX  = 


(3^)[(K-f)6rs  deRk 


rs 


+ 2G  de  ] 

rs 


(K  " T>  {W~  + 2G<W- 

rs  rs 


) - 


3F  3F 


3a  . p 
rs  3e 

rs 


(B-ll) 


(B-12) 


This  may  be  rewritten  as 

<;9F 

dX  = 


7T- — ) [Kde„„6  + 2G  de  ] 

3a  KK  rs  rs 


r s 


K( 


3a 


rs  rs 


3F 


3 3a 


pq 


5 <5 

pq 


r s 


2 

) - 


3F  3F 


(B-l 


30  3,-P 

rs  3e 

r s 


in  which  de  = de  - — d£i,i,6_0  are  the  deviatoric  strains, 
rs  rs  3 KK  rs 

The  incremental  stress  strain  relations  during  plastic  de- 
formation may  be  obtained  by  substituting  Eq.  (B-13)  into 
Eq.  (B-ll).  The  yield  surface  F(a^  > e?^ ) = 0 in  the  models 
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a.  . 


o. . 

13 


8ji  9oi3  2 vi ;avr2® 


6.  . + 


3ji  « avr: 


3F 

3ep 

rs 


3cP 
9F  v 


3F 


3eF 

v r s 


3e 


p r s 


Introducing  Eqs.  (B-18,  19)  into  Eqs.  (B-ll,  13)  gives 


3 F 3 r < a 

da..  = K5..(dekk  - 3dX  ~)  + 2G(de..  - dX — ~ ^ 

xj  H kk  3Ji  13  aVjJL  2H 


3F 

Tj”  Kde 
1 


kk 


+ G 


3 F 


r s 


dX  = 


3"V~J "^2 


9K( 


3F  ' 

3J1) 


+ G( 


3F 


3\T7 


-)  -3-^- 

' ’ 3J, 


de 


rs 


3F 

3ep 


For  the  case  of  uniaxial  strain  in  the  E-direction, 
following  relations  hold 

e = e = £ = e =e  =0 

xx  yy  xy  yz  zx 

a =0  = a 

xx  yy  r 

Denoting  the  axial  stress  and  strain  by  0 and  e , 


") 

2 


the 


respectively,  one  can  write 


dGkk 

= de 

z 

(B-24) 

de 

zz 

2 , 

= r de 
3 z 

(B-25) 

de 

XX 

* de  = - 4 de 

yy  3 z 

(B-26) 

de 

xy 

= de  = de  =0 

yz  zx 

(B-27) 

J1  = 

o +20 
z r 

(B-28) 

s = 

zz 

4(o  - a > 

3 z r 

(B-29) 

s = 

XX 

S = - 4(o  - 0 ) 

yy  3 z r 

(B-30) 

Sxy 

S = S =0 

yz  zx 

(B-31) 

[7~  - 

... 1 i a n i _ "VT"  | „ | 

(B-32) 

/J2  ' 

YJ  ' z r'  2 ^ z 

Substitution  of  these  equations  into  Eq.  (B-20,  21)  gives 
(3K  M + G — *’z  ) dG 

dA  1 9^j2  VT  |sj  2 

9K(^)2  + G(-^7"-33f  -2L 

9J1  ^fj2  9J1  3eP 


(B-33) 


and 


ds  = 2G  (—  de  - dX  s ) 

z 3 z 5 T z 

9J2 


(B-34) 


dJ.,  = 3K(de  - 3dX  •— — ) 
1 z 9J1 


Furthermore , 


d£ 


p - dX  r?~  «...  3dX  3F 


3°ij  13 


3J, 


(B-35) 


(B-36) 
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Equations  (B-33,  34,  35,  36)  may  be  integrated  numerically 
to  obtain  the  stress-strain  curve  of  the  model  in  uniaxial 
strain  provided  that  Eq.  (B-33)  can  be  evaluated  for  dX  . 
This  can  always  be  done  except  when 


3F  _ 3F 
8ji  ’ 


0 


(B-37) 


which  can  happen  for  the  model  of  Section  IV  only  for  F - F^ 
if  and  only  if 

S =J=£P=L=X=0  (B-38) 

2 1 V 

In  this  case  one  may  employ  an  asymptotic  approach  in  which 
all  quantities  which  vanish  are  expanded  about  the  origin  in 
stress  space 


J1 

= ji 

t 

dJl 

= ji 

dt 

(B-39) 

e 

= i 

t 

de 

= e 

dt 

(B-40) 

z 

z 

2 

z 

X 

= A 

t 

9 

dX 

= X 

dt 

(B-41) 

s 

= s 

t 

ds 

= s 

dt 

(3-42) 

2 

z 

z 

z 

ep 

♦ p 
= e„ 

t 

dep 

= ep 

dt 

(B-43) 

V 

V 

V 

V 

where  t is  a parameter  which  is 

to  approach  zero. 

Then 

3F9 

« • 

337  ■ 2(Ji  ' L>  ‘ 

2(Jl  - L) t 

(B-44) 

9f2  2 n~ 

3'/i;  v 2 

l*»l 

(B-45) 
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Substituting  Eqs.  (B--39  to  45)  into  Eqs.  (B-33  to  36)  gives 


Xt 


[6K.(J.,  - L)  + 2G.R2  s ]e 
x 1 i x z z 


36^^ 


L)  + 3G.R*  s2  -6  (J 
X X z 1 


L) 


3F 

3ev 


1 

t 


• O • • O • 

s = 2G . (—  £ - XR.  s t) 

z x 3 z x z 


J,  - 3K.[e  - 6X  (j,  - L) t ] 

i x z 1 

e£  = 6X  t (J,  - L) 


in  which  K.  , G.  and  R.  are  the  bulk  modulus,  shear  modulus 
xxx 

and  cap  shape  factor  under  conditions  (B-38) . 


Consider  Eq.  (16)  of  Section  IV, 

no 

i - X = R[A  - Ce  ] 


By  means  of  Eqs.  (22,  24)  it  is  clear  that  l > X , which 
implies  that  i > 0 . Therefore,  by  Eq.  (16),  L is 
instantaneously  equal  to  zero,  even  though  H is  not  zero. 
For  L = 0 


3F, 


3F 


3eP 


_2 

3x 


-S2L-  -2  (X  - L)-d* 


de 


dep 

v 


-2Xt 


d X 
dep 


= -2t( 


-±h  ep  - 

dep  v 
v 


-12Xt2  j (— ) 
1 dep 

v 


Then  Eq.  (B-46)  becomes 


Xt 


(6K , J . + 2G.RT  s ) £ 
i x x x z'  z 


. 2 4*2  . 2 d X • 

36K.J/  + 3G  Rj  s + 72  J.  (-2— ) Xt 
x i i i z 1 d£p' 

V 


(B-46) 

(B-47) 

(B-48) 

(B-49) 


(B-50) 


(B-51) 

(B-52) 
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s = ~r  » J = ~r  » A = * q ■ 

-z  ^z  def. 


(B-53) 


Then  Eqs.  (B-47,  48,  52)  become 


4G . „ 

S = —■  - 2G  A RT  s 
3 i l 


(B-54) 


J « 3K.  - 18  K.  A J 
l x 


(B- 55) 


6K.J  + 2G.RTS 

x 1 l 

36K.J2  + 3G.R*  s2  + 72  J2  A Q2 

ill 


(B-56) 


Equations  (B-54,  55,  56)  may  be  solved  for  s,  J and  A 
After  some  algebra  one  finds 


(M  _ 

. '■G.  3K . ' 

a = — — i i 

6 (4  - R25) 

3K. 

j = i 

1 + 18  K . A 

i 


(B-57) 


(B-58) 


S = E,J 


(B-59 ) 


in  which  5 = s/J  is  the  slope  of  the  stress  path  and 
satisfies  the  equation 


/r#  - iib  77  = (4  ' V1  + t R2i  *2 

i i de^ 


(B-60) 


The  initial  modulus  M for  uniaxial  strain  is  given  by 


• • 1 • 
s + i Ji 
M.  - 7^  = -Z-^-  3—1 
1 e,  e. 


s + k J = (€  + 


3(5  + §)  K.(4  - R2  5)  (5  + |)(4  - R25) 

4 ' * + 3K‘^  ' ^ <£-£->« 

3.  1 


(B-61) 
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Equation  (B-61)  is  useful  for  fitting  the  model  in  the  low 
stress  region  where  the  seismic  behavior  of  the  material 
may  be  important. 


